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I. INTRODUCTION 

The study of electron transport through a molecular 
system attached to metallic electrodes is regarded as 
one of the most promising research fields in nanoscale 
technology and physics^. With the discovery of ad- 
vanced molecular scale measurement methodologies like 
scanning electro-chemical microscopy (SESM), scanning 
tunneling microscopy (STM), atomic force microscopy 
(AFM), etc., it is now possible to measure current flow 
through single molecules or cluster of molecules sand- 
wiched between two electrodes^. The proposed idea of 
designing molecule-based diode by Aviram and Ratner— 
in 1974 first illustrates the possibility of using molecules 
as active components of a device. Since then several 
ab-initio as well as model calculations have been done 
to investigate molecular transport theoretically^"—. But 
experimental realizations took a little longer time to get 
feasible. In a pioneering experiment, Reed et at inves- 
tigated current-voltage {I-V) characteristics in a single 
benzene molecule coupled to metallic electrodes via thiol 
groupsi^. Various other experiments using molecules 
have also been reported in literature exploring many in- 
teresting features e.g., ballistic transport, quantized con- 
ductance, negative differential resistance (NDR), gate 
controlled transistor operation, memory effects, bistable 
switching, conformational switching to name a few. Al- 
though a lot of theoretical^"— and experimental22ri2S 
studies have been made so far using different molecules, 
yet several problems are to be solved for further develop- 
ment of molecular electronics. 

Most of the works associated with electronic conduc- 
tion through molecular bridge systems are mainly con- 
cerned on the overall conduction properties. But a very 
few works are available where attention has been paid to 
the current distribution within the molecule itself hav- 
ing single or multiple loop substructures^^"— . Recently 
some interesting works have been done by Nitzan et al. 
and few other groups where possible quantum interfer- 
ence effects have been explored on current distribution 
through such molecular geometries due to the existence 
of multiple pathways, yielding the possibilities of voltage 
driven circular currents^"—. The appearance of circular 



currents in loop geometries have already been reported 
in other contexts several years ago. This is commonly 
known as persistent currents in mesoscopic conducting 
rings where the current is induced by means of Aharonov- 
Bohm flux (j> threaded by the ring^"— . The reason be- 
hind this current in isolated loop geometries is quite dif- 
ferent from the previous one where current is driven by 
an applied bias voltage. It has been verified that the cir- 
cular currents appear in molecular rings, driven by ap- 
plied bias voltage, produce considerable magnetic fields 
at the centers of these rings. This phenomenon is some- 
what interesting and can be used in different aspects in 
the study of molecular transport. For example, in pres- 
ence of a local spin at the ring center one can regulate 
spin dependent transport through the molecular wire by 
tuning the orientation of that local spin, and, also the 
behavior of spin inelastic currents can be explained. In a 
recent work Galperin et al^ have proposed some results 
towards this direction. One can also utilize this circular 
current generated magnetic field in other way to control 
spin dependent transport through a molecular wire with- 
out changing the orientation of the local spin. In that 
case we can change the strength of the magnetic field 
by some mechanisms. To test it, biphenyl molecule may 
the best example where two benzene rings are connected 
by a single C-C bond. It has been examined^ that in 
the case of a biphenyl molecule, electronic conductance 
changes significantly with the relative twist angle among 
two benzene rings. For the planar conformation conduc- 
tance becomes maximum, while it gets decreased with in- 
creasing the twist angle. This phenomenon motivates us 
to describe conformation-dependent circular currents in 
a biphenyl molecule coupled to two metallic electrodes. 
We use a simple tight-binding (TB) framework to de- 
scribe the model quantum system and evaluate all the 
results through Green's function formalism. We believe 
that our present analysis will certainly provide some im- 
portant information that can be used to design molecular 
spintronic devices in near future. 

The structure of the paper is as follows. In section 
II, we describe the molecular model and theoretical for- 
mulation for the calculations. The essential results are 
presented in Section III which contains (a) transmission 
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probability as a function of injecting electron energy and 
junction current through the molecular wire as a function 
of applied bias voltage for different twist angles, and (b) 
conformation-dependent circular currents in two benzene 
rings and associated magnetic fields at the ring centers. 
Finally, in section IV, we summarize our main results and 
discuss their possible implications for further study. 



II. MOLECULAR MODEL AND THEORETICAL 
FORMULATION 

A. Tight-binding model 

Figured] gives a schematic illustration of the molecular 
wire, where a biphenyl molecule is coupled to two semi- 
infinite one-dimensional (ID) metallic electrodes, com- 
monly known as source and drain. Our analysis for the 
present work is based on non-interacting electron picture, 
and; within this framework, TB model is extremely suit- 
able for analyzing electron transport through a molecu- 
lar bridge system. The single particle Hamiltonian which 
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FIG. 1: (Color online). Schematic diagram of a biphenyl 
molecule attached to two electrodes, namely, source and 
drain. The magenta arrow describes the relative twist among 
the molecular rings. 

captures the molecule and side-attached electrodes gets 
the form: 



H = H. 



mol 



ele 



(1) 



The first term Hmoi corresponds to the Hamiltonian 
of the biphenyl molecule sandwiched between two elec- 
trodes. Under nearest-neighbor hopping approximation, 
the TB Hamiltonian of the molecule composed of 12 
[N = 12) atomic sites reads, 

i i 

■j j 
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where the indices i and j are used for the left and right 
molecular rings, respectively, e denotes the on-site en- 
ergy of an electron at i-(j-)th site and v describes the 
isotropic nearest-neighbor coupling between the molecu- 
lar sites, cj(cj) and Ci{cj) are the creation and annihi- 
lation operators, respectively, of an electron at the 
)th site. The last term in the right hand side of Eq. [5] 



describes the coupling between two molecular rings. In 
terms of the relative twist angle 9 among these two rings, 
the coupling strength Vi^r gets the form: 114.7 = vcos^ 9. 

The second and third terms of Eq. [1] describe the 
Hamiltonians for the ID semi-infinite electrodes (source 
and drain) and molccule-to-electrode coupling. In Wan- 
nier basis representation they are expressed as follows. 



Hele — Hs + Hd 



dn + h.c. 



a=S,D k " 



(3) 



and. 



.mol 

= Ts [cldo + h.c] + TD [c\dN+i + h.c] . (4) 

Here, eq and tg correspond to the site energy and nearest- 
neighbor hopping integral in the electrodes. d\^ and d„ 
are the creation and annihilation operators, respectively, 
of an electron at the site n of the electrodes. The coupling 
strength between the source and the molecule is rg, while 
it is td between the molecule and the drain. The source 
and drain are coupled to the molecule through p-th and 
g-th atomic sites, respectively, those are variable. 



B. Two-terminal transmission probability and 
junction current 

To obtain transmission probability of an electron from 
source to drain electrode through the molecule, we use 
Green's function formalism. Within the regime of coher- 
ent transport and in the absence of Coulomb interaction 
this technique is well applied. 

The single particle Green's function operator repre- 
senting the entire system for an electron with energy E 
is defined as, 



G= [E-H + iT])'^ 



(5) 



where, 77 —> 0'*'. Following the matrix forms of H and 
G, the problem of finding G in the full Hilbert space 
of H can be mapped exactly to a Green's function Q 
corresponding to an effective Hamiltonian in the reduced 
Hilbert space of the molecule itself and we have^, 



mol 



(6) 



Here, Us and So are the contact self-energies intro- 
duced to incorporate the effect of coupling of the molecule 
to the source and drain, respectively. In terms of this 
effective Green's function Q, two-terminal transmission 
probability T through the molecular wire can be written 
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where, Ts and To are the coupUng matrices, and, Q'^ 
and Q"" are the retarded and advanced Green's functions, 
respectively. 

With the knowledge of the transmission probability we 
compute overall junction current (It) as a function of 
bias voltage (V) using the standard formalism based on 
quantum scattering theory. 



function formalism. At absolute zero temperature the 
current flowing from site i to j (j = i ± 1) is given by 
the expression. 



lij — 



J .hj{E)dE, 



(11) 



It{V) 



2e 



T[fs{E)- fomdE. 



(8) 



Here, fs and fo are the Fermi functions of the source 
and drain, respectively. At absolute zero temperature the 
above equation boils down to the following expression. 



It{V) = ^ f T{E)dE, 



(9) 



where, Ep is the equilibrium Fermi energy. 



C. Circular current and associated magnetic field 

In order to calculate circular current in molecular rings 
of the biphenyl molecule let us first concentrate on the 
current distribution in a simple loop geometry illustrated 
in Fig. [2] A net current It flows between two electrodes 




FIG. 2: (Color online). Schematic view of current distribution 
through a ring geometry coupled to two electrodes. The filled 
black circles correspond to the positions of the atomic sites. 

through a quantum ring, where Ii and I2 are the currents 
flowing through upper and lower arms of the ring, respec- 
tively. We assign positive sign to the current propagating 
in the counter-clockwise direction. With this current dis- 
tribution, we define circular current of the ring as^S, 



(/iLi +/2L2) 



(10) 



where, Li and L2 are the lengths of the upper and lower 
arms of the ring, respectively, and L ~ Li + L2. Thus, in 
order to compute /j,, following the above relation, we need 
to know the currents in different branches of the loop 
geometry. We evaluate these currents by using Green's 



where, Jij is the current density. In terms of the corre- 
lation function it can be written as^^. 



(12) 



where, Q"' = Q^TsG'^- This correlation function is eval- 
uated by setting the occupation function of the source to 
unity and that of the drain to zero. 

Finally, we determine the local magnetic field at any 
point r inside the ring, associated with circular current 
Ic, from the Biot-Savart's Law^S, 
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dr X {r — r') 
' |(7^-r^)|3 



(13) 



where, fiQ is the magnetic constant and r' is the position 
vector of an infinitesimal bond current element lijdr'. 
Using the above expressions we evaluate circular currents 
and associated magnetic fields in two molecular rings of 
the biphenyl molecule. 

Throughout this work, we assume that the entire volt- 
age is dropped across the molecule-to-electrode interfaces 
and this assumption is good enough for molecules of 
smaller size. We also restrict ourselves at absolute zero 
temperature and choose the units where c = e = h = 1. 



III. NUMERICAL RESULTS AND DISCUSSION 

In this section we present numerical results computed 
for transmission probability, overall junction current and 
circular currents in a biphenyl molecule under conven- 
tional biased conditions. Throughout our analysis we 
set the on-site energies in the molecule as well as in the 
source and drain electrodes to zero, e = eg = 0. The 
nearest- neighbor coupling strength in the electrodes (to) 
is fixed at 2eV, while in the molecule (v) it is set at leV. 
The coupling strengths of the molecule to the source and 
drain electrodes, characterized by the parameters T5 and 
r£), are also set at leV. We fix the equilibrium Fermi en- 
ergy Ep aX zero and measure the energy scale in unit of 

V. 



A. Transmission probability and junction current 

In Fig. 13] we show the variation of transmission prob- 
ability T as a function of injecting electron energy E 
when the source and drain arc coupled to the biphenyl 
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molecule at the sites 1 {p = 1) and 10 {q = 10), re- 
spectively. The red curve corresponds to 6* = 0, while 
the green and blue lines arc associated with 6 ~ tt/S 
and 7r/2, respectively. From the spectrum it is observed 
that in the planar conformation (6 = 0) transmission 
probability exhibits resonant peaks (red curve) for some 
particular energies associated with the energy eigenval- 
ues of the molecule, and, at these resonances it (T) al- 
most reaches to unity. The situation is somewhat in- 
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FIG. 3: (Color online). Transmission probability as a function 
of energy for the biphenyl molecule when the electrodes are 
connected at the molecular sites 1 and 10, as shown in Fig.[T] 
The red, green and blue lines correspond to 6 = 0, n/3 and 
7r/2, respectively. 

teresting when we twist one benzene ring with respect 
to the other. It shows that the magnitude of the res- 
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FIG. 4: (Color online). Total transport current as a func- 
tion of bias voltage when the electrodes are connected to the 
molecule following the same configuration as prescribed in 
Fig. El 

onant peaks gets reduced compared to the planar con- 
formation with increasing the twist angle 6 and for the 
particular case when 9 becomes tt/2 i.e., two molecular 
rings are perpendicular to each other, transmission prob- 
ability exactly drops to zero (blue curve) for the entire 
energy band spectrum. With the increase of the twist 
angle between these two molecular rings, the degree of 
TT-conjugation between them decreases, which results the 
reduction of the junction conductance since the electronic 
transfer rate through the molecule scales as the square of 
the TT-overlap^. Thus, rotating one benzene molecule 



relative to the other one can regulate electronic trans- 
mission through the biphenyl molecule, which leads to 
the possibility of getting a switching action using this 
molecule. 

The sharpness of the resonant peaks in T-E spectrum 
strongly depends on the molecular coupling strength to 
side-attached electrodes, and, it controls electron trans- 
fer through the bridge system significantly. In the limit 
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FIG. 5: (Color online). Total transport current for a spe- 
cific voltage bias as a function of twist angle for the biphenyl 
molecule. The source and drain are coupled to the molecular 
sites 1 and 10, respectively. 

of weak-coupling i.e., ts{td) << v, sharp resonant 
peaks will be observed in transmission spectrum, whereas 
widths of these peaks get broadened in the limit of strong- 
coupling, ts{t£)) ~ V. The broadening of transmission 



















i 


i 








i 


I) 




-3 


-2 


-1 









2 


3 



E 



FIG. 6: (Color online). Transmission probability as a function 
of energy for the same parameter values used in Fig. [S] when 
the source and drain are coupled to the molecular sites 6 (p = 
6) and 9 {q — 9), respectively. 

peaks with the enhancement in coupling strength is quan- 
tified by the imaginary parts of the self-energy matrices 
Y.S and E^i which arc incorporated in the transmittance 
expression via the coupling matrices Fs and Td- This 
coupling effect on electron transport has already been 
explored elaborately in literature^"—, and therefore, in 
the present work we do not illustrate this issue further. 

The essential features of electron transport will be 
more transparent from our current-voltage characteris- 
tics. The overall junction current It through the molcc- 
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ular wire is computed by integrating over the transmis- 
sion curve, Eq. 1^1 following Landauer-Biittiker formal- 
ism. In Fig. S] we display the variation of junction cur- 
rent It as a function of applied bias voltage V for the 
biphenyl molecule where the electrodes are connected at 
the molecular sites 1 and 10, same as in Fig. [1] The 
current varies quite continuously with the voltage bias. 
Depending on the molecular coupling to the source and 
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FIG. 7: (Color online). Total transport current as a function 
of bias voltage when the electrodes are coupled to the molecule 
following the same configuration as mentioned in Fig. [5] 

drain electrodes, the current exhibits continuous-like or 
step-like behavior. Therefore, for a fixed voltage bias 
one can regulate the current amplitude by tuning molec- 
ular coupling strength, and, this phenomenon provides 
an interesting behavior in designing molecular electronic 
devices. From the spectrum (Fig. 2]) it is observed that 
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FIG. 8: (Color online). Internal current distribution in the 
molecule for its planar conformation when the bias voltage is 
fixed at 4V. The source and drain are attached to the atomic 
sites 1 and 10, respectively. 

the junction current decreases with increasing the relative 
twist angle, following the T-E characteristics. In addi- 
tion to this behavior, it is also important to note that the 
threshold bias voltage Vth of electron conduction strongly 
depends on the twist angle 0, which is clearly observed 
by comparing the red and green curves in Fig. |4l 

In order to explore the dependence of electron con- 
duction through the biphenyl molecule for any arbitrary 
angle of twist, in Fig. [5] we present the variation of 
junction current as a function of relative twist angle 9 
for some typical bias voltages. The spectrum reveals 
that for the planar conformation total current ampli- 
tude becomes maximum and it gradually decreases with 
the relative twist angle and eventually drops to zero 
when 6 = 7r/2. Thus, at 9 = 7r/2 no electron con- 



duction will take place through this molecular bridge 
system, while for other choices of 9 electron can trans- 
mit from the source to drain electrode, which promotes 
a conformation-dependent switching action using the 
biphenyl molecule. 

A significant change in the transmission spectrum is 
noticed when the electrodes are coupled to the molecule 
in such a way that the upper and lower arms of the molec- 
ular rings have unequal lengths. In Fig. [5] we present 
the results for such a particular configuration where the 
source and drain are attached to the molecular sites 6 
and 9, respectively. The red, green and blue curves cor- 
respond to the results for the same parameter values cho- 
sen in Fig. [3l From the transmission curves (red and 
green) we notice that for a wide energy range around 
iJ = 0, electron conduction does not take place, and also 




FIG. 9: (Color online). Internal current distribution in the 
molecule for its planar conformation when the electrodes are 
coupled to the atomic sites 6 and 9, for the same bias voltage 
taken in Fig. [S] The directions of the corresponding magnetic 
fields at the ring centers are illustrated by the encircled cross 
and dot representing downward (into page) and upward (out 
of page) directions, respectively. 

the widths of the transmission peaks get reduced enor- 
mously compared to the symmetric configuration where 
upper and lower arms of the molecular rings are identical 
(Fig. [31), even though the molecule-to-electrode coupling 
strength is unchanged. This is solely due to the effect of 
quantum interference among the electronic waves pass- 
ing through different arms of two molecular rings. In- 
troducing more asymmetry among these two arms one 
can regulate electron conduction through the bridge sys- 
tem in a tunable way. This behavior is nicely refiected in 
the current- voltage characteristics. The results are pre- 
sented in Fig. [T] It is clearly observed that for a partic- 
ular bias voltage the current amplitude decreases signifi- 
cantly compared to the symmetric configuration. Fig. 2] 
The symmetry breaking among the molecular arms also 
tunes the threshold voltage Vth of electron conduction 
across the molecular wire for a particular twist angle 0, 
as noticed by comparing Figs. |4] and [71 

B. Circular current and magnetic field 

Now we focus our attention on the behavior of circular 
currents and associated magnetic fields at the ring centers 
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of the molecule. These factors are highly sensitive to the 
molecule itself as well as the molecule-to-electrode inter- 
face geometry. To address these issues, we start with the 
current distribution within the molecule when it is cou- 
pled to the electrodes in such a way that the upper and 
lower arms of the molecular rings have identical lengths. 
The current distribution is shown in Fig. [SJ where the 
blue and green arrows indicate the bond currents in up- 
per and lower arms of the rings, respectively. The ar- 
row sizes represent the magnitudes of bond currents and 
they are computed when the bias voltage is fixed at 4V. 
Here It is the net junction current, shown by the red 
arrow, which is distributed among different branches at 
the junction point. Due to the geometrical symmetry 
reasons, the magnitudes of the bond currents in upper 
and lower arms of the two rings arc exactly identical and 
since they are propagating in opposite directions, no net 
circulating current will appear which results vanishing 
magnetic fields at the ring centers. 

In order to establish circular currents in these two 
molecular rings, we attach the electrodes asymmetrically 
such that the upper and lower arms of the rings have un- 
equal lengths. The internal current distribution for such 
a particular configuration is illustrated in Fig. [9l where 




currents and the associated magnetic fields at the ring 
centers as a function of voltage bias for the planar con- 
formation of the biphenyl molecule. From the spectrum 
we notice that in one voltage regime no circulating cur- 
rent appears, while other voltage regimes finite circular 
currents arc available and the sign of these currents also 
changes depending on the voltage region. The associated 
magnetic fields also follow the same behavior, as illus- 
trated in the spectrum. This phenomenon can be ex- 
plained as follows. The circular current in a loop geome- 
try is associated with energy eigenstates those can be de- 
scribed by current carrying states with the current flow- 
ing in opposite directions. Now, for a finite bias voltage 
whenever one of these resonant states lies in the Fermi 
window, associated with the applied voltage bias and the 
nature of voltage drops along the molecular wire, we will 
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FIG. 11: (Color online). Circular currents and corresponding 
magnetic fields at the ring centers for a specific voltage bias 
{y — 1.5V) as a function of twist angle when the molecule is 
sandwiched between electrodes according to the configuration 
taken in Fig. [9] The red and black lines correspond to the 
identical meaning as in Fig. 1101 



FIG. 10: (Color online). Circular currents and associated 
magnetic fields at the ring centers of the biphenyl molecule as 
a function of bias voltage when the electrodes are coupled to 
the molecule following the configuration prescribed in Fig. [51 
The red and black curves describe the results for the left and 
right rings, respectively. The twist angle 9 is fixed at zero. 

we set the same bias voltage as taken in Fig. [51 In this 
situation the bond currents get unequal magnitudes, and 
accordingly, circular currents are established in the two 
molecular rings. The magnetic fields at the ring cen- 
ters associated to these circular currents are also shown. 
In the left ring, the magnetic field is directed along the 
downward direction, represented by encircled cross, while 
in the other ring it is directed along the upward direction, 
illustrated by encircled dot. We calculate the magnetic 
fields using Biot-Savart's Law, Eq. [T31 and scale them in 
unit of Gho/AttR, where R is the perpendicular distance 
from the center to any arm of the molecular ring and the 
factor 6 appears due to the existence of six bonds in each 
benzene ring. 

Figure [TUl demonstrates the magnitudes of the circular 



get the corresponding circular current. When more than 
one resonant states come within this Fermi window, all 
of them contribute to the current which provide a net 
signal, and, in the particular case when they mutually 
cancel each other, the net signal becomes zero. The sign 
of this net circular current or direction of associated mag- 
netic field depends on which resonant states dominate the 
others, which again certainly depends on the applied bias 
voltage. 

Finally, we focus on the variation of circular currents 
and corresponding magnetic fields at the ring centers as 
a function of relative twist angle 0, when the voltage bias 
is kept constant. The results are shown in Fig.[TTl where 
two different colored curves represent the similar meaning 
as in Fig. 1101 Quite interestingly we see that the magni- 
tudes of the magnetic fields at the centers of two molecu- 
lar rings decreases monotonically with relative twist an- 
gle 9, and for large enough 9 they eventually reduce to 
zero. Thus, for a particular bias voltage one can tune the 
strength of magnetic field established at the ring center 
due to this circular current simply by rotating one ben- 
zene molecule relative to the other, and hence, by placing 
a local spin or a magnetic ion at the ring center, which 
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will interact to the magnetic field, spin selective trans- 
mission can be achieved through this molecular system. 
This conformation-dependent spin selective transmission 
will be investigated in a recent forthcoming paper. 

IV. CONCLUSION 

In conclusion, we have investigated in detail the 
conformation-dependent two-terminal electron transport 
through a biphenyl molecule within a simple tight- 
binding framework using Green's function formalism. 
Two principal results have been obtained and analyzed. 
First, the dependence of molecular twist on electronic 
transmission probability and the overall junction current 
have been discussed. Our results lead to a possibility of 
getting the conformation-dependent switching action us- 
ing this molecule. Second, we have investigated the vari- 
ation of circular currents and associated magnetic fields 
developed at the ring centers as a function of the relative 
twist angle. Tuning this angle, one can tailor the strength 



of magnetic field at the ring centers, and, we believe that 
the present analysis may provide the possibilities to de- 
sign molecular spintronic devices using organic molecules 
with loop substructures. 

Throughout our work we have ignored the inter- and 
intra-site Coulomb interactions as well as the effect of the 
electrodes, which we plan to consider in our future works. 
Another important assumption is the zero temperature 
approximation. Though all the results presented in this 
communication are worked out at absolute zero temper- 
ature limit, the results should remain valid even at finite 
temperatures (~ 300 K) since the broadening of the en- 
ergy levels of the biphenyl molecule due to its coupling 
with the metal electrodes is much higher than that of the 
thermal broadening^. 
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